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Abstract 


We propose a model for addressing the superfluidity of two different Fermi species confined in 
a bilayer geometry of square optical lattices. The fermions are assumed to be molecules with 
interlayer s—wave interactions, whose dipole moments are oriented perpendicularly to the layers. 
Using functional integral techniques we investigate the BCS-like state induced in the bilayer at 
finite temperatures. In particular, we determine the critical temperature as a function of the 
coupling strength between molecules in different layers and of the interlayer spacing. By means of 
Ginzburg-Landau theory we calculate the superfluid density. We also study the dimerized BEG 
phase through the Berezinskii-Kosterlitz-Thouless transition, where the effective mass leads to 
identify the crossover from BGS to BEG regimes. The possibility of tuning the effective mass as a 
direct consequence of the lattice confinement, allows us to suggest a range of values of the interlayer 
spacing, which would enable observing this superfluidity within current experimental conditions. 
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Introduction.- Although the underlying mechanism for the occurrence of high Tc super¬ 
conductivity has not been fully elucidated, it has been established that it must include 
two essential features, long range interactions among its constituents and transport in par¬ 
allel layers [l, 2|. At present, experiments carried out with ultracold molecules offer the 
possibility of emulating such conditions as a result of their precise control achieved in the 
laboratory Indeed, successful realizations of ultracold dipolar gases of molecules, with 

magnetic or electric nature, either with permanent or induced dipole moment have already 
been performed. Besides this interest, it has been proposed that efficient quantum informa¬ 


t.on p.oce=.„g a can be achieved wbh ayaten. w.th long cange natecact.ona n_t 
conhned in a 2D geometry. Similarly, the predicted Wigner crystalline phase [10| also shares 
those distinctive attributes. 


Here, we address on the study of superfluidity in dipolar Fermi molecules conhned in 2D 
optical lattices. For this purpose we consider the system described below and work within 
the mean-held perturbative analysis using the integral functional formalism to determine 
the critical transition temperature. The superhuid transition is analyzed in the Ginzburg- 
Landau (GL) scheme where the superhuid density is calculated. We found that assisted by 
the 2D lattice conhnement the attractive interactions between molecules in diherent layers 


induces interlayer paring 
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13| of both types, BGS and bound molecular BEG states. 


Model- Our system consists of a gas of dipolar Fermi molecules placed in two parallel 
optical lattices in 2D. The conhguration of the optical lattices in both layers has the same 
structure, a square lattice with constant lattice of size a. See Fig. [H In the presence of an 
electric held perpendicular to the layers, the dipoles are aligned along the same direction 
and, consequently, their interactions within the same layer are suppressed by intermolecular 
repulsion |3|. In contrast, dipoles in diherent layers attract each other at short range, 
while repelling each other at large distances. Thus, the interaction among dipoles has a 
specihc form resulting from such a conhguration. As schematically shown in Fig. [H such an 
array can be mapped into a system of fermions in two diherent hyperhne spin states, with 
interactions within an ehective 2D environment. Fermions moving in up and down layers 
can be associated to two diherent hyperhne spin states labeled with 'I' and respectively. 
The interaction potential between fermions in diherent layers (or diherent hyperhne spin 
states) has the form 111 
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FIG. 1. (Color online) Schematic representation of the dipolar Fermi gas. Polar molecules in the 
up (down) layer can be mapped into the specie labeled with t (),) when the gas is described in a 
2D layer. 



( 1 ) 


being r the in-plane separation between fermions, d the dipole moment and A the interlayer 
spacing. The Hamiltonian that describes such a system is 



H 


0(=A,B 


( 2 ) 


where HqIv) is the ideal term that includes the kinetic energy and the external potential 
created by the 2D optical lattice V)ati(T) = Vq (sin^(a;7r/a) -|- sin^(r/7r/a)). The held operators 
^a(f) = = J2k^ki^')^k satisfy the usual commutation relations for 


fermions. The in-plane energy dispersion of the ideal Fermi gas within the tight binding 
approximation is = —2f(cos k^a -|- cos kyo) being t the hopping strength among near¬ 
est neighbors. It is worth mentioning that in real experiments, besides the optical lattice 
conhnement, the atoms typically move under the inhuence of a harmonic conhnement. How¬ 
ever, since the frequency uj of this potential (due to the magnetic trapping) is so small, the 
curvature imposed on the optical lattice can be neglected; that is, tvjj t. Therefore, the 
Hamiltonian of the two component Fermi gas in the momentum representation adopts the 
form 



k 


k,k',q 
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where Vdipik — k') is the Fourier transform of Vdipif — f'). In the last equation we have 
introdnced the chemical potential to acconnt for the conservation of the total molecule 
number, and have already considered the lattice geometry. We note that hi = Nx x Ny, 
being Nx and Ny the nnmber of sites along x and y directions respectively. 

To stndy the snperflnidity in the model, we use functional integral techniqnes. The action 
for the ideal Hamiltonian term is, 

So = ^ <p\k) ■ [ihwni - (eg - fi)a^] ■ 4>{k), (4) 

k 

where we use the abbreviations, k = {k,Un), 0 = {(l>A{k),(j)%{k))^ with 0 q, Grassmann 
fermionic nnmbers, I is the identity in 2D, and the fermionic Matsnbara freqnencies are 
hun = i^n + l)7r//9, with (3 = l/ksT. a a (a = x,y,z) are the Pauli matrices. To write 
the action associated to the interaction energy term, we introdnce the auxiliary Hubbard- 
Stratonovich transformation in terms of the bosonic held A 

(A(fc)) = - >• ( 5 ) 

k' 

that is, the action of snch held creates a pair of particles with the same momentnm k in 
diherent layers A and B. In terms of the helds 0 and A the action is 


(e^ - + A(k, g)ax) -^(k, g)—^ ^ A{g, k)V^A{k-k')A{g, k'), 

k^q k,k',q 

( 6 ) 

being <F(/c, g) = {(l>A,k, 4’B,q-kY'■ Since this action has a qnadratic form in the helds $ one 
ends with the following expression for the ehective action 


S.„[A] = ^lVln (Cy) 

k^q 


^ k'), 

k,k',q 


(7) 


with = (ihcunl — {ej: — jj,)az + A{k, g)ax). From this eqnation and within hrst order 
pertnrbation theory, the eqnation for the critical temperatnre, that signals the transition to 
the snperhnid state, is 


^o{k) = ^ Vdipik - k 
k' 


y Ao(fc')tanh {l3^p/2) 


2e 


( 8 ) 


k' 
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where = y + Ao(A;)2. In this work we have restricted to the case g = 0 and 

to an s-like gap. That is, assuming the momentum of the center of mass of the two-body 
collision term being zero, as in the original BCS scheme. 

Critical Temperature. - The physics of the system is completely determined by the dipolar 
interaction strength ad = meff(P/h^ and the dimensionless parameters A = A/a and y = 
a^/A, with rUe// = /2ta? as the effective mass. In principle, within our scheme, interlayer 


BCS pairs are formed for arbitrary values of these parameters 


14l |. On the other side. 


consideration of the two-molecule collisions is in order, since, as in the homogeneous case 
[n|. the two-molecule interlayer potential Vdip leads to different regimes of scattering as 
a function of the molecular dipole moment and the interlayer spacing. Namely, bound 
molecular states or BCS pairing can occur. Below we discuss the region of the parameters 
at which bound states can be formed. 
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FIG. 2. (Color online) Energy gap vs thermal energy (Eq. ([8])) for A = 0.5. Calculations were 
done for square lattices of = Ny = 120 and half filling occupation. 


To determine the critical temperature as a function of the parameters described above, 
we assume half-hlling occupation. For all of our numerical calculations we work within the 
first Brillouin zone (—^ < k^^ky < and use lattices of size AA = Ny = 120, which 
produce the same quantitative results than bigger lattices. The chosen values for molecular 
dipole moments and masses correspond to those used in recent experiments on ultracold 
polar molecules 1^. In Fig. (|2]), we can observe the critical temperature at which the 
superfluid state occurs, within the BCS regime, for different values of the coupling strength 
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X and fixed A = 0.5. Further, in Fig. |3]we plot the BCS critical temperature Tc versus x 
for different values of A, in order to illustrate the dependence of the critical temperature on 
those parameters. 
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FIG. 3. (Color online) Critical temperature as a function of the dimensionless interaction strength 
X for A = 0.5, 0.75,1.0. 


Ginzburg-Landau Free-Energy.- We now expand the effective action, given by Eq. ([7]), in 
terms of the new bosonic held A around the phase transition, and relate it to the Ginzburg- 
Landau Free-Energy Q 

S,„lA,A-] = ft/Jfi(A). (9) 

To do this, we write G~^^{k,q) = Go^(fc) -|- 'E^{k,q), where GQ^(fc) = ihwrJ. (Jz{ek — /a) 
andE^{k,q) = axA(k,q), thus G~^^{k,q) = GQ^{k){I-\-Go{k)T.^(k,q)). From the expansion 
lnG“^(/c,g) = lnGg^(/c) -|- X]m=i m we verify that odd powers of GqS do 

not contribute to the trace in S'e//[A], while even powers are explicitly given by 


JoSa)"’" = I- 




22 ) 


1 2m ’ 


[{hwn)^ + (e^-/r)2]' 

being and G 22 = — — As usual within the GL theory and 

well justihed near the pairing transition temperature (A small), we truncate the expansion 
up to quartic order. This procedure leads us to determine the coefficients of the free energy 
F = a(r)A^ -|- 5(T)A^. As expected, the coefficient a{T) bears the same information as 
that given by the gap in Eq. ([8]), while b{T) conveys the information on the stability of the 
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equilibrium state. We emphasize that within the mean field approach, the whole information 
about the dipole interaction is contained in the critical temperature Tc. The coefficient b{T) 
is found to be, 

+ -^——^l- (10) 


w=^5:1-/5^ 


4 (efc-/r)3 


It is important to emphasize that, while the expression for b{T) appears to be independent 
of the 2D lattice confinement and of the two-molecule interaction potential, their influence 
is exhibited through the critical temperature Tc, which it has been shown above to, in turn, 
strongly depend on the dipolar potential, namely on y and A. As expected, the coefficient 
b(Tc) is always positive. 

Superfluid density.- To determine the superfiuid density we also work within the per¬ 


turbative scheme, neglecting quantum fluctuations, and following 


Landau free-energy theory for second order phase transitions 
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he functional Ginzburg- 


2l|. That is, we expand 
the free energy in terms of the order parameter A, incorporating a phase A = Aoe*®'^% 
where 0 = ]^)- This procedure allows us to separate the density into the normal and 

superfiuid components. In order to make explicit the dependence on the phase parameter 0 
and using the symmetry associated to the conservation of total molecule number, a gauge 
transformation can be performed on the operators that create and annihilate particles at 
the lattice sites in both layers: (A|, Afl, {Bj, Efl (.B|e*®'^% 5*6“®'^*). 

By performing this transformation in Hamiltonian given by Eq. Dll. having written it in 
its lattice representation 22|], one ends with the following expression for the Hamiltonian 

Dll: 


He = H„ 


E(4“£+4A) E (® 


a=x,y 


der 0? dh 

+ 


'dkc 


2 dkl 


+ Hj, 


( 11 ) 


that is, associated to the phase Q ■ Ri, there appears a contribution to the kinetic energy. 


Thus, in this scheme, the superfiuid density component can be calculated as 




Pa,a' = lim 


1 Fp 




0 


0—>-0 A^f 0rv0n'^ 


Nt dOrydOr 


-,a= {x.t/}, 


( 12 ) 


being Fq and Fq the free energies of the normal and superfiuid phases respectively. The 
determination of the superfiuid density pa,a'i using Eq. flT^ . can be straightforwardly ac¬ 
complished by noting that in the effective action, see Eq. o. an extra term appears asso¬ 
ciated to the phase of the gap. Thus, one can write, in an explicit way: G^^(/c) —)■ G^^^k) 
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and consider pertnrbations np to second order. That is, keeping qnadratic terms of 0 in 
Ge^(fc) = GQ^Q{k) (I — G 0 =o(A:)S©(/c)). The hnal expression for pa,a' is 

1 dh. 


paa' — 


Ga2 


^ (n(fc) 


dkndkni 


_ Y(k)^^ 

' ’ak„ak„, 


(13) 


where n{k) is the momentnm distribntion and Y(k) is the Yoshida distribntion dehned as 
Y(fc)=/3sech^(/3^^/2)/4, with = \/(e^ — pY + Since the off-diagonal terms of the 
snperflnid tensor are small p = {p^x + Pyy)/‘^- FignrelH shows the snperflnid density fraction 



FIG. 4. (Color online) Superfluid density for different values of the dimensionless interaction 
strength y and A = 0.5. 


as a function of the temperature for three different values of the interaction strength. We 
denote by po the superfluid density fraction at T = 0. As it is well known the superfluid 
density fraction at zero temperature is different from unity in the lattice 211 . 

Berezinskii-Kosterlitz-Thouless temperature.- When the interaction y is strong enough, 
the predicted critical temperature becomes an artifact of the approximation since the BCS 
approach fails. To determine the relationship of the critical temperature with the coupling 
interaction, the presence of bounded pairs should be included. This procedure is delineated 


in Ref. 


23| and gives rise to a clear identihcation of the Berezinskii-Kosterlitz-Thouless 


(BKT) transition temperature, as that at which a fraction of bound pairs dissociate. Thus, 
to estimate the region at which our treatment is valid, we compare it with the mentioned 
procedure. This yields the following result. 

Tip 


ksTsKT — 


8m, 


eff 


(14) 
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In Fig. [5] we plot the critical temperature as a function of the dimensionless interaction 



FIG. 5. (Color online) BCS critical temperature and the Berezinskii-Kosterlitz-Thouless Tqkt 
transition temperature as a function of the interaction strength x, for A = 0.5. The inset correspond 
to the binding energy of the two-molecule system in the presence of the lattice. 


strength x for both, our BCS scheme and the incorporation of bounded pairs, see Eq. 
O- We observe how in the weakly interacting regime the temperature Tbkt is close to 
TbcSj while they separate in the strong interaction region y > 0.3. Beyond this value, the 


critical temperature Tbkt corresponds to a phase in which dimerization dominates 


2^. It 


is important to note that while in the homogeneous system the mass is hxed, here, the 
presence of the lattice determines the effective mass, as a tunable parameter, that yields 
the scale at which the transition from BCS phase to dimerized BEC phase occurs. The 
inset of Fig. |5] shows the pair binding energy using a variational wave function for the two- 
molecule problem, ^(r, 7 ) = with 7 a variational parameter, within the presence of the 
2D lattices. 

The system of dipolar Fermi molecules in a 2D square lattice here proposed can be 
realized as an application of the recent experiments with ultracold molecules with anisotropic 
interactions. In particular, those with with polar molecules ^°K®^Rb conhned in optical 
lattices, of a lattice constant of 520 nm approximately 15|. For molecules such as ^°K®^Rb 
and ®Li^^^Cs the dipole strength d can be varied with an external electric held from 0.2 
to 1.5 D. Those experiments allow us to suggest typical interlayer spacings of A between 
200 — 500 nm, well within current experimental limits. 
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In conclusion, we have studied superfluidity in polar Fermi molecules placed in 2D square 
optical lattices bilayers. We have used a functional integral formulation to study the finite 
temperature BCS state as a function of the coupling strength y and the dimensionless pa¬ 
rameter that measures the scaling between the interlayer separation and the lattice size. 
To determine the superfluid density we employed the Ginzburg Landau free energy scheme 
by exploiting the symmetry associated to the isotropy in the lattice and the total molecule 
number conservation. Given the current experimental context, our study constitutes a strat¬ 
egy for handling the critical temperature of the superfluid state by changing both the dipole 
strength d and the interlayer spacing among the 2D square optical lattices. 
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